* Replicate Figure 1 *
use "KS_state_level.dta"
twoway (scatter dpctgop casratepermil, mlabel(stab) msymbol(none) xlabel(10 15 20 25 30 35 40 45) ylabel(-.15 "-15%" -.10 "-10%" -.05 "-5%" 0 "0%" .05 "5%" .10 "10%"))  (lfit dpctgop casratepermil, ytitle("Change in GOP Vote Share, 2012 to 2016" " ")  xtitle(" " "State Casualties per Million Residents") graphregion(color(white)) leg(off))

* Replicate Table 1 *
use "KS_county_level.dta"
replace medianfaminc = medianfaminc/1000
reg dpctgop2 casrateper10k 
outreg2 using table1, word label dec(3) replace
reg dpctgop2 casrateper10k pctcollege medianfaminc pctwhite pctblack pctlatino pctrural pctvets
outreg2 using table1, word label dec(3) append
